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Abstract —Sharing radar spectrum with communication sys¬ 
tems is an emerging area of research. Deploying commercial 
wireless communication services in radar bands give wireless 
operators the much needed additional spectrum to meet the grow¬ 
ing bandwidth demands. However, to enable spectrum sharing 
between these two fundamentally different systems interference 
concerns must be addressed. In order to assess interference 
concerns we design a three-dimensional (3D) channel model 
between radar and cellular base stations (BSs) in which the radar 
uses a two-dimensional (2D) antenna array and the BS uses a 
one-dimensional (ID) antenna array. We formulate a line-of-sight 
(LoS) channel and then propose an algorithm that mitigates 
radar interference to BSs. We extend the previously proposed 
null space projection algorithm for 2D channels to 3D channels 
and show that effective nulls can be placed by utilizing both the 
azimuth and elevation angle information of BSs. This results in 
effective interference mitigation. In addition we show that the 3D 
channel model allows us to accurately classify the size of radar’s 
search space when null space projection algorithm is used for 
interference mitigation. 

Index Terms —MIMO radar, null space projection, spectrum 
sharing, channel modeling, full dimensional MIMO, 3D beam¬ 
forming, 3D channel. 

1. Introduction 

Today’s wireless communication operators urgently require 
more spectrum to increase capacity and meet consumer de¬ 
mands. A classic practice to free up more usable spectrum 
for wireless systems has been to relocate incumbents to some 
other band. However, not only this process is time consuming 
it is also very expensive for incumbents to relocate. An 
emerging way forward, to avoid the above mentioned issues, is 
to share spectrum across government agencies and commercial 
services. In the United States, the Federal Communications 
Commission (FCC) is considering to allow commercial wire¬ 
less operation in the 3.5 GHz radar band Q. This initiative 
has the potential to leverage new methods of spectrum sharing 
among heterogeneous systems such as small cells, radars, 
and satellite systems. However, a fundamental limitation for 
spectrum sharing is interference caused to coexisting systems. 
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Novel interference mitigation methods are required for coex¬ 
istence of these fundamentally different systems. This letter 
addresses this issue by modeling a 3D channel between radar 
and communication systems and proposes an algorithm that 
mitigates radar interference to cellular BSs. 

A fundamental area of research in wireless communications 
is channel modeling as it allows performance evaluation of 
transmission techniques. Traditionally, wireless channel mod¬ 
els have been designed to capture the azimuth i.e. 2D chan¬ 
nels. A comprehensive introduction to 2D wireless channel 
modeling, propagation modeling, and statistical description of 
channels can be found in iQ and references there in. Recently, 
there has been a surge in modeling wireless channels that 
capture both the azimuth and elevation directions known as 
3D channels. 3D channels allow very precise beamforming 
for users on ground and in buildings. As a result the wireless 
research community is heavily involved in designing and 
standardizing 3D channel models that can reap the benefits 
of both the azimuth and elevation beamforming O . However, 
the efforts so far have been limited to model 3D channels 
between wireless communication systems and to the best 
of our knowledge no work exists on 3D channel modeling 
between wireless communication and radar systems. A 3D 
channel model can enable radar systems to place accurate nulls 
in azimuth and elevation angles of BSs in order to mitigate 
radar interference. Therefore, in order to study performance 
evaluation and interference mitigation techniques in a spec¬ 
trum sharing scenario, the problem of 3D channel modeling 
is of prime importance. 

In this letter, we formulate a 3D MIMO channel model 
between MIMO radar and MIMO cellular system. We consider 
a 2D antenna array at the radar and a ID array at the BS. Using 
our proposed channel model we demonstrate the efficacy of a 
novel null space projection algorithm in 2D which mitigates 
radar interference. The rest of this letter is organized as 
follows. Section [III briefly presents MIMO radar architecture 
and spectrum sharing scenario. Section jllll models 3D MIMO 
channel between seaborne radar and cellular system. Section 
HVl proposes a novel interference mitigation algorithm for 
3D channels. Section El discuss the volume illuminated by 
spectrum sharing MIMO radars that are subject to interference 
mitigation. Section [^discusses the simulation results. Section 
IVIII concludes the letter. 


11. Spectral-Coexistence Models 


In this section, we briefly introduce the fundamentals of 
MIMO radar and discuss our radar-communication system 
spectrum sharing scenario. 

A. MIMO Radar 

We consider a shipborne mono-static MIMO radar with 
M = Mh X My colocated antennas in a uniform rectangular 
array (URA), where is the number of horizontal antenna 
elements in a given row and My is the number of vertical 
antenna elements in a given column. We assume that the inter¬ 
element distance between adjacent antenna elements on any 
given column and row is dg . The target is assumed to be a point 
target which has azimuth and elevation angles represented by 
0 and (j), respectively. 

Let X be the transmitted waveform vector which simul¬ 
taneously illuminates the target area of interest and nulls 
interference to BSs present in a 2D spatial sector. 

B. Spectral Coexistence Scenario 

We consider a spectrum sharing scenario in which a cellular 
system is deployed in the radar band. Each BS of the cellular 
system is equipped with N transmit and receive antennas in a 
uniform linear array (ULA). Each BS antenna is at a speciflc 
azimuth and elevation angle with respect to radar. We assume a 
Erequency Division Duplex (EDD) cellular network in which 
BSs are operating in the radar band and UEs are operating 
at non-radar frequencies. Therefore, we focus on interference 
caused by the radar operation to the base stations only. We 
consider a single cell or a BS with % users, without loss of 
generality, that receives the following signal in the presence 
of radar 

X 

r = ^GiSi + Hx + n (1) 

i=l 

where is the channel gain between the BS and the i* user, 
H is the channel gain between the BS and the radar, and x is 
the interfering signal from the MIMO radar. This interference 
from radar can be mitigated by projecting radar waveform onto 
null space of interference channel between radar and BS PI- 
lEl. Previous work considers null space projection for 2D 
channels O which is extended in this paper to 3D channels. 

III. 3D Channel Modeling 

In this letter, we formulate a 3D channel model between 
radar and communication system. We introduce a 2D uniform 
rectangular array (URA) antenna array at the radar and a ID 
uniform linear array (ULA) at the BS. We assume M = M^ x 
My antenna elements, where M^ antenna elements exist in 
the horizontal/row direction and My antenna elements exist in 
the vertical/column direction with equal element spacing dg. 
Each element of the 2D antenna array can be labeled by s, 
where s = (/, k) = {k — l)My + /, for /c = 1, 2,..., Mh and 
/ = 1,2,...,M^. 

We proceed by deriving 3D channel model for line-of- 
sight (LoS) case. In the absence of any beamforming, the 
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Eig. 1. Spectrum sharing scenario between radar and cellular BSs. The 
portion containing BSs is nulled-out to protect it from the radar’s interference. 


3D channel characteristics are dependent on element radiation 
pattern, antenna spacing, and propagation of each ray path. In 
the LoS case, for the two adjacent vertical antenna elements 
s and s + My from the same column, the channel coefficients 
between receive-transmit antenna pairs {u, s) and (fx, s + M^), 
i.e., Hu s and Hu,s-\-Mv, have the following property 

H^,s+M^ « cos(^) (2) 


i.e. the difference between the two adjacent elements can be 
captured by the phase shift due to propagation. In equation 
©, 0 is the angle of departure in elevation domain. 

Eor a BS with N antenna elements in a uniform linear array 
(ULA), the N x MMIMO channel is partitioned into M^ sub¬ 
matrices given by 


H^ = 



iik=2 




(3) 


where each has size N x My and is given as 


H^ = 


^2,(/c-l)M^+2 

^Ar,(/c-l)M^+2 




for /c = 1,2,... ,M/,. 

Using equation (O we can write the above channel matrix 
as 






tik=Mh 


^ih{O)0lik=l (4) 


where ( 8 ) is the Kronecker product and 


a/,(6>) = [1 


^-j27r/X{Mh-l)ds cos(0)j 


is the steering vector in azimuth domain. Using the same 
arguments we can write the partitioned channel model in 
elevation domain as 




H;=i 


H;=m„ 


(5) 












and each H; has size N x M„. Using equation Q we can 
write above as 




H;=1 Hi=2 


H;=m„ 


a„((/))®H,=i (6) 


where 


av{(f>)=[l e--^ 27 r/Ad, cos(<^) 


^-j 27 r/X(M^-l)ds cos((/))j 


is the steering vector in elevation domain. 


IV. 2D Null Space Projection 

In this section, we introduce NSP in both azimuth and 
elevation domains by considering the 3D channel model de¬ 
signed in the previous section. We define projection algorithm 
which projects radar signal onto null space of 3D interference 
channel, i.e., in elevation and azimuth H/^. Assuming, the 
MIMO radar has channel state information, we can perform 
singular value decomposition (SVD) in both domains to find 
the null space and then construct a projector matrix. We 
proceed by first finding SVD of H^, i.e.. 


2; 



Eig. 2. Search space of spectrum sharing MIMO radar. The part containing 
cellular BSs is not part of the radar’s illuminated area due to the nulls placed to 
protect BSs from its interference. That is why the blocked-out part (blue part) 
of the figure is an illustration of the nulls placed in that particular azimuth- 
elevation sector. 



(7) 

Now, let us define 


= diag((j7;^i, c^^!,2 5 • • • 5 ^v,p) 

(8) 

where p = mm{N,My) and d-yp > ay^2 > 

> (Jy^q > 

o'v,q-\-i = <^v,q -\-2 = ''' = = 0 are the singular values of 

H^. Next, we define 


^v,2i • • • 5 ^v,My 

) (9) 

where 


A JO, for u<q, 

11, for u > q. 

(10) 

Using above definitions we can now define 
matrix in elevation, i.e.. 

our projection 

P„ 4 . 

(11) 

Using similar arguments we can define projection matrix in 

azimuth, i.e.. 



(12) 


where Heiv and HAzm afe waveform covariance matrices in 
elevation and azimuth, respectively. The 3D beampattern of 
NSP signal is given as 

G{0,(j)) = {(/)) Knsp (18) 

V. Radar Search Space 

Radar search space is normally specified by a search solid 
angle O in steradians. If we define the radar search volume 
extent for both azimuth and elevation, in degrees, as 0 and 0, 
the search volume can be given as EH 

" = ( 57 ^ (stetadians)^ ( 19 ) 

Let OnuII denote the sector, azimuth and elevation angles, of 
cellular BSs that needs to be nulled to protect BSs from radar 
interference. Then, due to NSP the search volume, shown in 
Figure 2, can be given as 

Onsp = ^(57 296)2* (steradians) ( 20 ) 


Using the above defined projection matrices we can write 
covariance matrices of projected waveform in elevation as 

(13) 

and azimuth as 

= P^R^,„PT (14) 

The covariance matrix of NSP transmitted waveform is given 


by 

^NSP = '^v,NSP + ^h,NSP 

(15) 

where 

^v,NSP = ^Azm^h{0)ay {(I))'ReIv 

(16) 


^h,NSP = KAzmM0)^vW^ErJ' 

(17) 


VI. Numerical Examples 

In this section, we provide numerical examples of 3D chan¬ 
nel modeling between radar and communication system and 
its applications using null space projection. For the simulation 
setup we assume that the radar has a 40 x 25-element URA 
and the BS has a 10-element ULA. We assume a target is 
located at 0 = 0° and (j) = 50° and communication systems 
are located in a sector of azimuth angles 0 = —45° to —40° 
and elevation angles 0 = 5° to 15°. We intend to null this 
sector to mitigate radar interference to this area. We provide 
two figures of the same example with different views to better 
understand the achieved results. 

In Figures [3] and [H the impact of 2D NSP is studied on 
the radar beampattern. It can be noted that the target can 


















Fig. 3. A 3D view of interference mitigation at BSs. Target is located at 
0 = 0° and = 50° and BSs are located in a sector of azimuth angles 
6 = —45° to —40° and 0 = 5° to 15° and we intend to null this sector to 
mitigate radar interference to this area. 

be accurately localized when it is far from the nulled areas. 
Moreover, the power received at the target is 60 dB which 
is significantly higher than the power received at the BSs 
which is less than -40 dB in this example. Thus, the proposed 
approach satisfies both the radar mission objective of target 
detection and simultaneously mitigates interference to BSs. 
Moreover, only the areas where BSs are present are nulled. 

In Figure [5l we show the volume illuminated by the spec¬ 
trum sharing radar as a function of its distance from nulled 
BSs. We assume the search-able space is 180° in azimuth, 
i.e., from -90° to 90° and 110° in elevation, i.e, from -20° 
to 90°. Note that for close-up range the difference in search 
volume with NSP and without NSP is significant, i.e., for ships 
docket at the harbor or close to shore line will experience 
some degradation in search volume due to NSP. However, if 
the radar is more than 2000 meters away from BSs the loss 
in search space is minimal. 

VII. Conclusion 

Spectrum sharing among heterogeneous systems is a 
promising solution to solve the spectrum congestion problem. 
In this letter, we explored a spectrum sharing scenario between 
seaborne radar and on-shore cellular system. We formulated a 
3D channel model between radar and communication system. 
We assumed a 2D antenna array at the radar and a ID antenna 
array at the BS. In addition, we introduced a new interference 
mitigation algorithm that projected radar waveform onto null 
space of 3D channel to mitigate radar interference to cellular 
BSs. The efficacy of the proposed algorithm was shown using 
simulation results. Search space of spectrum sharing MIMO 
radars using NSP was also explored and a relation to quantify 
the size of search space was provided. It was shown that the 
reduction in search space was directly related to the areas 
nulled to protect BSs. 



Azimuth 0 

Fig. 4. An image of interference mitigation at BSs. The intensity/magnitude 
of radar signal can be observed at the target and BSs. Note the difference 
between the power received at target vs. the power received at BSs. 
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Fig. 5. Search volume for radar located far from BSs. Search-able space is 
180° in azimuth, i.e., from -90° to 90° and 110° in elevation, i.e, from -20° 
to 90°. So the searchable volume in percentage at a distance of: 500 m is 
95.9%, 2 km is 99%, and at 8 km is 99.8%. At longer distances the search 
volume with and without NSP is very close. 
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